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Abstract 

This article presents a unified approach to simultaneously compute the Jacobians 
of several singular matrix transformations in the real, complex, quaternion and 
octonion cases. Formally, these Jacobians are obtained for real normed division 
algebras with respect to the Hausdorff measure. 

1 Introduction 

A fundamental tool in statistical theory and in particular in distribution theory is the com- 
putation of Jacobians of matrix transformations. Many such Jacobians were first found i n 
the real case and subsequently f or comp l ex cas e s; among many o t hers, see Woodind ( 1956 ): 
IJamesI (Il964i Sect ions 4 and 8); iKhatril (|l965l ). iMuirheadl (|l982h ; iRatnaraiah et all (j2005l) 
and lMathail (|l997t) for a detailed review of this topic. Most of these Jacobians were obtained 
with respect to the Lebesgue measure. However, several recent articles have examined these 
Jacobian for the singular random matrix case, i.e., they study densities and then calculate 
the Jacobian s of the matrix transformations with res p ect to the Hausdorff measure, se e 
Khatril (Il968l). lUhlig I (Il994l). iDi'az-Garcia et all (Il997l). Ipiaz-Garci'a and Gutierrezl (Il997l). 



Srivastaval(l2003l)TDiaz-Garcia and Gutierrez- Jaime^ "(l2006h. l Diaz-Garcfa and Gonzalez-Farias! 



(|2005al) . lDiaz-Garcia and Gonzalez-Farias I (l2005bl). lip et al\fa007\ ) and lDiaz-Garcial^OOTl ) 



among others. In addition, IRatnaraiah and Villancourti (120051) st udied some of these lat- 
ter Jacobians in the singular complex case and iLi and Xue (2010) considered the singular 
quaternion case. 

Using results obtained from abstract algebra, it is possible to demonstrate a unified 
means of addressing the computation of Jacobians in the singular and nonsingular real, 
complex, quaternion and octonion cases, simultaneously. This approach has been used for 
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Edelman and Raol (|2005h 



som e time in r andom matrix theory in the nonsingular case, see 
and lForresteH (|2009h . 

For the sake of completeness, in the present work we consider the case of octonions, but 
it should be noted that many results for the octonion case can only be c onjectured, because 
there remain many unresolved theor etical p roblem s in this respect, see Drav and Manoguel 
(1999). furthermore, as stated by iBaed (|2002h . the relevance of the octonion case for 
understanding the real world has yet to be clarified. 

The rest of this paper is structured as follows: Section [5] provides some definitions and 
notation on real normed division algebras, showing some ideas about the explicit form of 
the Hasusdorff measure. The main results on the computation of Jacobians for singular 
matrix transformations are presented in Section [3] It is emphasised that all these results 
are obtained for real normed division algebras. 



2 Notation and real normed division algebras 



Let us introduce some notation and useful results. A detailed discussion of real normed 



division algebras may be found in iBaea (|2002l ). For convenience, we shall introduce some 
notation, although in general we adhere to standard forms. 

For the purposes of this study, a vector space is always a finite-dimensional module 
over the field of real numbers. An algebra g" is a vector space that is equipped with a 
bilinear map m : $ x $ — > 5 termed multiplication and a nonzero element 1 £ J termed the 
unit such that m(l,a) = m(a, 1) = 1. As usual, we abbreviate m(a, b) — ab as ab. We do 
not assume g associative. Given an algebra, we freely think of real numbers as elements of 
this algebra via the map w t— > uil. 

An algebra g is a division algebra if given a,b g # with ab = 0, then either a = or 
6 = 0. Equivalently, ^ is a division algebra if the operation of left and right multiplications 
by any nonzero element is invertible. A normed division algebra is an algebra # that is 
also a normed vector space with ||a6|| = ||a||||6||. This implies that J is a division algebra 
and that ||1|| = 1. 

There are exactly four normed division a lgebra s : real numbers (5ft), complex numbers 
(<£), quaternions (i}) and octonions (D), see iBaez (|2002h . We take into account that !K, 
£, f) and O are the only normed division algebras; moreover, they are the only alternative 
division algebras, and a l l divis ion algebras have a real dimension of 1,2,4 or 8, which is 
denoted by /3, see IBaez ( 20021 Theorems 1, 2 and 3). In o ther branches of mathem atics, 
parameters a = 2/(3 and t = /3/4 are used, see Tableland Edelman and Raol ( 2005 ) and 
Kabd (1 19841) . 



Table 1: Values of f3 = 2/ a and t = /3/4 parameters. 



p 


a 


t 


Normed divison algebra 


1 


2 


1/4 


real (5ft) 


2 


1 


1/2 


complex (£) 


4 


1/2 


1 


quaternionic (fj) 


8 


1/4 


2 


octonion (O) 



Let n be the linear space of all n x m matrices of rank m < n over $ with m distinct 
positive singular values, where J denotes a real finite- dimensional normed division algebra. 
Let g rnxm be the set of all n x m matrices over The dimension of g'™ xm over 5ft is j3mn. 

T 

Let A e g nxm 7 then A* = A denotes the usual conjugate transpose. 
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The set of matrices Hi G £™ xm such that H^Hi = I m is a manifold denoted V^„, 
and is termed the Stiefel manifold (Hi is also known as semi- orthogonal (0 = 1), semi- 
unitary (0 = 2), semi-sym vlectic (0 = 4) and semi- exceptional type (0 = 8) matrices, see 
Drav and Manoguel ( 19991) ). The dimension of „ over 3? is [/3mn — m(m — l)/3/2 — m]. 



In particular, m with dimension over 3?, [m(m + 1)0/2 — m], is the maximal compact 
subgroup il^(m) of £^ and consists of all matrices H G gmxro guc j 1 ^. nat _ j m _ 

Therefore, U^m) is the real orthogonal group 0(m) (0 = 1), the unitary group U(m) 
(0 = 2), compact symplectic group Sp(m) (0 = 4) or exceptional type matrices Oo(m) 
(0 = 8), for # = 3t, £, or £>, respectively. 

We denote by the real vector space of all S G 5 mxm such that S = S*. Let 
be the cone of positive definite matrices S G grmxm. th en is an open subset of 
6^. Over 3?, ©^ consist of symmetric matrices; over £, Hermitian matrices; over Sj, 
guaternionic Hermitian matrices (also termed self-dual matrices) and over O, octonionic 
Hermitian matrices. Generically, the elements of 6^, are termed Hermitian matrices, 
irrespective of the nature of The dimension of ©^ over 5ft is [m(m — 1)0 + 2] /2. 
Let ID^j be the diagonal subgroup of £^ m consisting of all D G g -mxm ; D = diag(di, . . . , d m ) 
and let Tu^ n be the semi-upper-triangular subgroup of £^ n consisting of all T G g -nxm j 
with > 0. 

Now, let £^f n (q) be the linear space of all n x m matrices of rank q < min(rt, m) with q 
distinct singular values and let &^f(q), the (0mq — 0q(q — l)/2) — q-dimensional manifold 
of rank q of positive semidefinite matrices S G g™ xm with q distinct positive eigenvalues. 

For any matrix X G J 1 "'™, dX denotes the matrix of diffe rentials (dxg), a nd we denote 
the measure or volume element as (dX) when X G g rmx ™ ) see iDimitriu ( 2002 ). 



If X G !fi nxm then (dX) (the Lebesgue measure in g nxm ) denotes the exterior product 
of the 0mn functionally independent variables 

n m ft 

(dX) = A A dxij where dXij = /\ dxy ^ . 

i=l j'=l fc=l 

If S G ©^ (or S G In+" n ) then (dS) (the Lebes gue measure in ©^ or in < Xu^^ n ) denotes 
the exterior product of the m(m + 1)0/2 functionally independent variables (or denotes the 
exterior product of the m(m — 1)0/2 + n functionally independent variables, if su G 3? for 
all i = 1, . . . , m) 

m p 

A A *C 

(dS) = { 4 - J 

f\d Sii f\ /\ds§\ if Sll G3?. 

i—l i<j k—1 

The context generally establishes the conditions on the elements of S, that is, if Sij £ 
G <£, G or G D. It is considered that 

(««) = A A*5 )s A*«A A*5 } - 

i<j fc— 1 i—l i<j fc— 1 

Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S G ^3^,, 
that is, S must be a non singular Hermitian matrix (Hermitian positive definite matrix). 
If A G 23^ then (dA) (the Lebesgue measure in 2)^,) denotes the exterior product of the 
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(3m functionally independent variables 



(dA) = A A dA 



i=i fe=i 

In addition, observe that, if X G £+^(5), we can write X as 



Xi 



Xn X12 

q X g qxm — q 

X21 Xn 

n—qxq n—qxra—q 



such that r(Xn) = g. T his is equivalent to the righ t product of the matrix X with a 
permutation matrix II, see Golub and Van Loan ( 19961 section 3.4.1, 1996), that is Xi = 
XII. Note that the exterior product of the elements from the differential matrix dX is not 
affected by the fact that we multiply X (right or left) by a perm utation matrix, that is, 
(dXj ) = (d(X Il)) = (dX), since n G UP(m), see iMuirheadl (Il982l Section 2.1, 1982) and 
Jamesl(ll954h . Then, without loss of generality, (dX) will be defined as the exterior product 



of the differentials dxij, such that Xij are mathematically independent. It is important to 
note that we will have (nq + mq — q 2 )f3 mathematically independent elements in the matrix 
X G C^f n (q), corresponding to the elements of Xn,Xi2 and X21. Explicitly, 



n q P q m 

(dx) = (dxu) a (dx 12 ) a (dx ai ) = A A A dx if AAA dx f 

i=lj=lfc=l i=l j=q+l k=l 



(1) 



Similarly, given S G &+?{q), we define (dS) as 



(dS) 



A A A ds 

1— l j=i fe=l 
q q m f3 

A^A A A ds 

2— 1 2—1 J=i+1 fe=l 



(k) 

ij J 



if Si, G 5ft. 



The context generally establishes the conditions on the elements of S, that is, if s,j G 5ft, 
G £, G S) or G £). It is considered that 

q m /3 q q m /3 

( dS ) = AA A<— A^*A A A d i fe) - 

i=lj=i/c=l i=l i—lj—i+lk—1 

Again, we should note that, for this case, the matrix S can be written as 



Sn S12 

qxq qxm—q 

S2I S22 

m—qxq m—qxm—q 



with 



r(Su) = q. 



such that the number of mathematically independent elements in S is j3mq — q(q + 1)/3/2 — q 
corresponding to the mathematically independent elements of S12 and Sn. Recall that 
Sn G in such a way that Sn has q(q — l)/3/2 + q, therefore, 

(dS) = (dS u ) A (dSu). 

Observe that an explicit form for (dX) and (dS) depends on the factorisation (base 
and coordinate set) employed to represent X or S, that is, they depend on the measure 
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factorisation of (dX) and (dS). For example, by using the nonsingular part of the decom- 
position in singular values and the nonsingular part of the spectral decomposition for X 
and S, then we can find an explicit form for (dX) and ( dS) (see Propositions 15711 and 15721 
respec tivel y), which are not unique, s ee Khatri ( 19681 ). IDiaz-Garcia et alH (|l997| ). lUhlig I 
( 19941 ) and lDiaz-Garcia and Gutierrez ( 1997| ). Alternatively, an explicit form for (dX) and 
(dS) can be found in terms of the QR and Cholesky decompositions, see Propositions 13.31 
and !3.4i respectively. 

A singular random matrix X in Z m ^ n {q) or &^f(q) does not have a density with respect 
to Lebesgue's measur e in g r " xm , but it does possess a de n sity o n a subsp ace M. C 5" xm ; see 
^ (l968h . lRao1 (I1973L p. 527). IDiaz-Garcia et al. I (Il997l ). lUhlig I (Il994h 



Khati 



and Cramer 



(1999, p. 297). Formally, X has a density with respect to Hausdorff's measure, which coin- 



cides w ith Leb esgue's measu r e, when the latter is define d on the subspace M; see|B illingslcy 



1986 . p. 2 47).IUhlig1 dl994l).lDiaz-Garcia et al. I dl997l ). lDiaz-Garcia and Gutierrez- Jaimez 
2005 ) and Diaz-Garcia and Gutierrez- JaimezT i 20091 ) 



The surface area or volume of the Stiefel manifold V£, „ is 



Vol«„) = / (H^dH 



/3/2] 



(2) 



where H = (H1IH2) = (hi, . . . , h m |h m +i, . . . , h n ) € il^(m). It can be proved that this 
differential form does not depend on the choice of the H2 matrix. When m = 1; vf n defines 
the unit sphere in g"™. This is, of course, an (n — 1)0- dimensional surface in g™. When 
m = n and denoting Hi by H, (H*dH) is termed the Haar measure on il^(m). Also, T^Ja] 
denotes the multivariate Gamma function for the space & m , and is defined by 



where etr(-) = exp(tr(-)) 
Gross and Richardd (19871) . 



/ ebr{-A}\A\ a -( m - 1 W 2 - 1 (dA) 

m 

n n(™-iw*Y[r[ a -(i-i)p/2], 

i=l 

denotes the determinant and Re(a) > (m — l)/3/2, see 



3 Jacobians 



We now consider several Jacobians of singular matrix transformations in terms of the B pa- 

I — I \ k I ' — 

ramet er . For a detailed disc u ssion of related issues see Uhlig (1994), Di'az-Garci'a and Guti errez 

(Il997l). IDiaz-Garcia et"aL~l dl997l), IRatnaraiah and Villancourtl ([20051 ) . lip et al\ (|2007l ). - 
Diaz-Garcia and Gutierrez- Jaimez (|2012l ) and lLi and Xud (|2oToh 



Pro positio ns 13.11 and 13.21 a nd Co r ollary 13.11 g eneralise the result s in IDiaz-Garcia et al 



(|l997l ). IRatnaraiah and Villancourtl (|2005l ) md Li and Xuel ( l20ic j) obtained for the real, 
complex and quaternion cases, respectively. 



Proposition 3.1 (Singular value decomposition, SVD). Let X 6 £^f n (q), such that 
X = ViDWJ be the nonsingular part of the SVD with, Vi G V^ n , Wi G 
D = diag(di, • • ■ , d q ) G D^, di > • • • > d q > 0. Then 



Vim and 



(dX) 



wild. 



■P(n+m-2q+l)-l 



JJ(d? - djf{dB) A (VJdVi) A (WJdWi) 



(3) 
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whe 



0, = 1; 
-j3q/2, = 2,4,1 



Proposition 3.2 (Spectral decomposition, SD). Let S 6 &^f(q). Then the nonsingular 
part of the spectral decomposition can be written as S = WiAWJ, where Wi £ Vq m and 
A = diag(Ai, . . . , A m ) € £>J, with Ai > • • • > A 9 > 0. TTiera 



(dS) = 2~% T f[ Af (m ^ 9) n(A, - Aj)"(dA) A (WJrfWi), 
where r is defined in Provosition \3. 1\ 

Corollary 3.1. Let X e £+f„(g), and S = X*X = WiAW^ £ 6+^g). Then 
(dX) = 2-«|A| /? ("- m+1 )/ 2 - 1 (dS) A (VJrfVi), 

with Vi G V£ ? . 

Proof. The proof is obtained immediately from Propositions 13.11 and 13.21 □ 

Let A 6 2>n^rn (l) then A + £ £+ ^ (g) denotes its Moore- Penrose inverse. 

Theorem 3.1 (Moore-Penrose inverse S = S*). Let S £ &£f(q) as in Proposition \3. 
Then ignoring the sign, i/V = S + £ (q) 



(4) 



(5) 



□ 



(dV)=HA 



/3(-2m+ g +l)-2 



(dS). 



(6) 



Proof. The proof follows by observing that if S = WiAWj is the nonsingular part of the 
spectral decomposition of S, then V = S + = WiA _1 Wj;. Therefore, from (0} 



(dS) = 2"% T fl Af (m ~ 9) JJ(Ai - A,y(dA) A (WJrfWi). 



And 



(dv) = 2-^Y[\: i 



-0{m-q)-2 



n(V 1 -A7 1 )^(W 1 dW 1 )A(dA), 



(7) 



taking into account that (ignoring the sign), 

dX 



(dA- 1 ) = /\ dAr 1 = A(-!)^ = n v 2 a dA< = n v 2 (^a). 



i=l 



(=1 



»=1 i=l 



(.=1 



Then, 



(W*dWi) A (dA) = 2%~ T 



-0{m-q)-2 



(dV). 



Finally (ignoring the sign), observe that, 



(Ai - A.,-)* 3 '-V A, ; 



n 



-0(9-1) 



The desired result is obtained by substituting © and © in ([7]). 



□ 



(8) 



(9) 
□ 
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Theorem 3.2 (Moore-Penrose inverse). Let X G £+^(<j<) as in Proposition \KT\ The 
ignoring the sign, ifY = X + 



(dY) = l[d: t 



-2/3(m+n-q) 



(dX). 



(10) 



j=l 



Proof. The proof is analogous to that given for Theorem 13.11 using ^ . □ □ 

Theorems l3.1l and l3.2l were obtained bv lDiaz-Garcia et al. I (1997) and Li and Xuel ( 2010( ) 
for the real and quaternion cases, res pectivel y. 

The n ext result was proposed by lUhligl (1994) in the real case as a conjecture. Sub- 
seque ntly, Di'az-Garci'a and Gutierrez (fl997 ) pro posed an indirect proof of this conjecture. 
Later iDiaz-Garcfa and Gutierrez- Jaimezl (|2009bf ) p rovided an a l ternat ive proof based on the 
exterior product, also in the real case. In 2010, iLi and Xud (l2010f) extended this result 
to the quaternion case, generalising the indirect proof stated in Dfaz-Garcia and Gutierrez! 
(jl997l ) . We now propose two alternative statements of these results for real normed division 
algebras. 

Theorem 3.3 (First Uhlig's conjecture via SVD). Let X,Y € &£f(n) such that X = 
B*YB, where B G £^ lm (m). In addition, consider the nonsingular part of the SD, X = 
GiAGJ and Y = HiAHJ, where A = diag(£i, . . . , S n ), A = diag(Ai, . . . , A„) <E T)\ and 
Gi.Hi G V? m . Then 



(dX) = |B|^"|G*B*Hi|' 3(m -"- 1 ) +2 ( ( iY) 
= |B^ n |H;BGi| fl ( m - n - 1 J +a (dY) 
= \Bf n \ A|^ (m_ri_1) / 2+1 |A| _ ' 3(m_ ™ _1)/2_1 (cZY), 



(11) 
(12) 

(13) 



with 



(dY) 



i—1 i<Cj 

Proof. Considering Propositions l3.1ll3.2l and lDiaz-Garcia and Gutierrez- Jaimezl (2012, Propo- 
sition 1) the proof is analogous to that given in iDiaz-Garcia and Gutierrez- Jaimezl (|2009bh 
for the real case. □ 

Proposition 3.3 ( QR decomposition, QRD). Let X G Z^ m (q), then there exists a matrix 
Hi G Vq n and an upper quasi-triangular matrix T G %u^ q such that X = HiT is the 
nonsingular part of the QR decomposition and 

1 



(dX) = Y[t^ n - l+1) - 1 (H* 1 dU 1 ) A (dT). 



(14) 



Proposition 3.4 (Cholesky's decomposition, CHD). Let S G 6+f(q). Then S = T*T, 
where T G c £u^ i q such that 

T=(TiT 2 ), 

with Ti G Tc/+^, q x q upper triangular matrix. Also, let X G S, m , n > with X = HiT (QR 
Decomposition) and S = X*X = T*T (Cholesky decomposition) such that 



S = 



T * Ti T * T2 



T^Ti 



T*T 2 



Sll S12 

qxq qxm — q 

S2I S22 

m—qxq m — qxm — q 



with 



Sn G^. 



Then 
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(dS) = 2 q ~[[t 



f)(m-i)+l 



(dT) 



(15) 



2!. ,4Zso, 



(dX) = 2-«|TJTi|' 3 ( n - m+1 ^ 2 - 1 (dS) A (HJdHi), 
= 2-«|S 11 | /3 ("- m+1 )/ 2 " 1 (dS) A (HJdHi. 



(16) 
(17) 



Propositions 13.31 and 13.41 were studied by Diaz-Garcia and Gonzalez-Farias ( 2005bl ) for 
the real singular case and bv lLi and Xuel (I2009T for the quaternion nonsingular case. 

We now study an alternative approach to Theorem 13.31 with respect to an alternative 
factorisation measure based on QR and the Cholesky decomposition. 

Theorem 3.4 (First Uhlig's conjecture via QRD). Let X, Y G 6+^(n), such that X = 
B*YB, with B G £+f m (m) fixed. Additionally, let X = T*T and Y = L*L with T,L G 
lo^, smc/i £/ia£ T = (Ti T2) and L = (Lj L2), where T\ and Li are nxn upper triangular 
matrices. Then 



(dX) 



TTT 



1 J-H 



m — n— 



n-l)/2+l| L * Li |-^( m -n-l)/2-l| B |^n( (f y^ 



(18) 



with 



(dY) = 2" JJ i 



/9(m— 



(dL). 



ZB. 



Proof. Let Z G £+^(n), such that Y = Z*Z. Then 

X = B*YB = B*Z*ZB = with * 

from Proposition ^. 4[ equation 1171 

(d$) = 2-"|TJT 1 |' 3 ("- m+1 ^ 2 - 1 (Q*dQ)(dX). 
In which * = QT with T G £tf+f n , Q G U p (n), and X = = T*T. Then 
(dX) = 2™|TtT 1 r /3 (™-" i+1 )/ 2 " 1 (Q*dQ)- 1 (d*). 



(19) 
(20) 
(21) 



Note that d<& = dZB, and then by Diaz-Garcia and Gutierrez- Jaimezl (|2012 , Proposition 1) 
we have that (d<&) = |B|^"(dZ), from which, substituting in ([2"Tj). we obtain 



(dX) = 2"|T*Ti|-' 3 ( n - m+1 )/ 2 - 1 (Q*dQr 1 |B|' 3n (dZ). 
Now Y = Z*Z = L*L with L G T+{*„; from Lemma EU 

(dZ) = 2-™|L*Li| /3 ("- m+1 >/ 2 - 1 (dY)(H*dH) 



(22) 



(23) 



where Z = HL, with L G c Zv' r ^ n and H G iiP jn). Moreove r, note that, due to the uniqueness 
of Haar's measure, (Q*dQ) =' (H*dH), see Ijamesl (|l954[ ). Thus, substituting (23) in (22), 
we obtain 

(dX) = |LtLi| /3( " _m+1)/2 ~ 1 |T*X 1 |"' 3 ( TI ^ m+1) / 2 - 1 |B|' 3 "(rfY) 
= |TtTi| /3(m - n+1)/2 |LJLi|' 3 - (m - n+1)/2 |B|' 3,l (dY). 



□ 



□ 
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The following result combines Theorems 13.11 and 13.31 which enables us to study the 
multivariate F and beta distributions (also termed matrix multivariate beta type II and 
type I, respectively) for real normed division algebras in a class of singularity. 

Theorem 3.5. Let X,Y G &tn( n ) such that x = B*Y+B, where B G £^ m (m). In 
addition, consider the nonsingular part of the SD, X = GiAG^ and Y = HiAH*, where 
Gi,Hi G V£ m and A = diag(<5i, . . . , 5 n ),A — diag(Ai,.. . , A„) G V\. Then 

(dX) = | B |' 3 "|A| /3 ( m -™- 1 )/ 2+1 |A|-' 3 ( 3m - n - 1 ^ 2 - 1 (dY), (24) 

with 

n n 

(dY) = 2- n n T Y[ Af [m - n) Y[(Xi - Xjf(dA) A (H*e£Hi). 

z=l i<j 

Proof. Let Z = Y+ and observe that Z = HiA _1 H^. Then by Theorem[S31 

(dX) = | B |' 3 "|A|' 3 ( m -"- 1 )/ 2 + 1 | A -i|-/»(«-»-D/2-i ( dz) . (25) 
Now by Theorem 13. II 

n 

(dZ) = Y[ \P { - 2m+n+1) - 2 {dY). (26) 

i=l 

The result follows by substituting JH into (gS]). □ □ 

Observe that this result corrects an err atum in the exponent of the determinant of A in 



Diaz-Garcia and Gutierrez- Jaimea (|2009bl . Theorema 2.2), obtained in the real case. 



Conclusions 

Note that the results presented here contain as special cases all versions of the results 
previously obtained in real, complex and quaternion and octonion cases, both for singular 
and nonsingular cases. For example from Theorem 13.41 

(rfX) = (|T 1 T 1 ||L 1 L 1 |- 1 ) /3(m_ ™~ 1)/2+1 \Bf n (dY), (27) 

Now, let X and Y be nonsingular matrices, that is, m = n, therefore Ti and Li are square 
nonsingular triangular matrices and 

|TJTi| = |X| = |B*YB| = |B*L^LiB| = |B| 2 |L^Li| 

then, 

\B\ 2 = \T* 1 T 1 \\L* 1 L 1 \- 1 . (28) 

The desired result 

(dX) = \B\^ m -^ +2 (dY), 
is obtained by sub stituting (|2"5)l into ([2"T]) , see Mathail ()l997l ) for real and complex cases and 



,mea by su 

Li and Xuei (|2010t ) for the quaternion case, among others. 



Finally, observe that the real dimension of real normed division algebr as can be expr essed 
as potentia of 2, /3 = 2™ for n — 0,1,2,3. Moreover, as observed by iKabd fll984h . the 



results obtained in this work can be extended to hypercomplex cases; that is, for complex, 
bicomplex, biquaternion and bioctonion (or sedenionic) algebras, which of course are not 
division algebras (except the complex algebra), but are Jordan algebras, together with all 
their isomorphic algebras. Note, too, that hypercomplex algebras are obtained by replacing 
the real numbers with complex numbers in the construction of real normed division algebras. 
Thus, the results for hypercomplex alge bras a r e obta ined by simply replacing j3 with 2/3 in 
our results (we reiterate, as reported by IKabd (|1984I )V 
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